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Abstract 

In this paper we complete the analysis begun by two of the authors in a 
previous work on the discrete quantum walk on the infinite line [J. Phys. A 
36:8775-8795 (2003); quant-ph/0303105 . We obtain uniformly convergent 
asymptotics for the "exponential decay" regions at the leading edges of 
the main peaks in the Schrodinger (or wave-mechanics) picture. This 
calculation required us to generalise the method of stationary phase and 
we describe this extension in some detail, including self-contained proofs 
of all the technical lemmas required. We also rigorously establish the 
exact Feynman equivalence between the path-integral and wave-mechanics 
representations for this system using some techniques from the theory of 
special functions. Taken together with the previous work, we can now 
prove every theorem by both routes. 
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1 Introduction 

The first authors to discuss the quantum walk were Aharonov, Davidovich and 
Zagury in 0] where they described a very simple realization in quantum optics. 
In this model a particle takes unit steps on the integers at each time step, 
starting at the origin. In j^H] Meyer proved that an additional spin-like degree 
of freedom was essential if the behaviour of the system was to be both unitary 
and non-trivial. Without this degree of freedom, the only way the evolution of 
the walk can avoid being purely ballistic is to relax the unitarity condition. This 
spin- like degree of freedom is sometimes called the chirality, or the coin, which 
is why this type of walk is sometimes called a "coined" walk This is in 

sharp contrast to the continuous time walk ^] ^1 which does not need a coin 
and which we will not discuss here. The chirality can take the values RIGHT 
and LEFT, or a coherent superposition of these. For a detailed introduction to 
quantum walks, we refer the reader to the review article in |21| . 

Meyer and subsequent authors |3 E] have considered two approaches to 
the discrete-time quantum walk, the path-integral approach of Feynman and 
the Schrodinger wave- mechanics approach, which reflect two complementary 
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Figure 1: The discrete quantum walk on the line. The probability distribution is 
shown for a walk that started at the origin with its coin in the state |L) after it has 
evolved for 100 steps. The distribution is oscillatory between the two peaks and decays 
exponentially outside that range. The peaks move away from the origin with speed 
t/y/2, and the width of the peaks also decreases with time. By contrast, the classical 
random walk has a Gaussian distribution, which spreads with velocity ~ \/t. 

ways of formulating quantum mechanics |15j . We refer to the paper by Ambai- 
nis, Bach, Nayak, Vishwanath and Watrous |S] for proper definitions and more 
references. Both approaches are discussed in the paper of Ambainis et al. The 
probability distribution for this walk is shown in Figure ^ after enough time 
has elapsed for the asymptotic behaviour to manifest. 

This paper began as a sequel to the work of Carteret, Ismail and Richmond 
concerning the one-dimensional quantum walk on the integers, and contains 
the completion of the analysis of this quantum walk for the remaining exponen- 
tial decay region in the Schrodinger picture. Wc believe the analysis presented 
in this paper to be interesting for three reasons. One is methodological; while 
analysing this system we encountered various links between a number of dif- 
ferent methods in combinatorics which do not seem to be widely known, and 
which may be of use to the quantum information community when analysing 
more complicated systems than the one discussed here. 

The second motivation is rather more abstract. It is one of the fundamen- 
tal principles of quantum mechanics that the wave-mechanics and path-integral 
representations of a system should produce exactly the same results. The quan- 
tum walk has been proposed as the quantum analogue of the classical random 
walk, in the hopes of ultimately defining a systematic procedure for "quantiz- 
ing" classical random walk algorithms |21l \G\ \'6'6\ . Quantizing classical systems 
is something that must be done with considerable care; the obvious approach 
isn't necessarily the correct one. While it is true that when such pathologies 
have been discovered in the past, they were found in much more exotic systems 
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than this one, it is as well to check. One way to perform such a check is to verify 
that the Feynman equivalence principle still holds between the path-integral and 
wave-mechanics representations. As it happens, the results obtained from the 
two approaches do not at first appear to be equivalent. In fact they are, though 
the proof of their equivalence is nontrivial, and will be given in Section|3| In the 
course of this analysis we uncovered a small, but potentially significant omission 
in previous analyses of this system, which we will describe below. 

Another reason for performing this check is to explore the little mystery 
left at the end of JUj. While the results from the path integral calculation 
made intuitive sense, the partial results from the wave-mechanics calculation 
for the exponential decay region were rather unexpected. Specifically, we found 
what appeared to be evanescent waves in this exponential decay region, which 
seemed to imply the presence of some kind of absorption mechanism, despite 
the fact that the definition of the model precludes any barrier or other source 
of dissipation in the system to cause these by the familiar mechanisms. So, we 
will check that the Feynman equivalence holds to verify that those evanescent 
waves are not simply some kind of mathematical artefact. 

We would also like to gain some insight into the physical interpretation of the 
mathematical behaviour of this part of the wave-function. A link between the 
behaviour of the quantum walk on the line and certain phenomena in quantum 
optics has been suggested previously by Knight, Roldan and Sipe in [2511241125) . 
and refined by Kendon and Sanders in 23 . This connection will be discussed 
in more detail in Section 0] 

We will also describe a potentially useful method for obtaining integral rep- 
resentations of orthogonal polynomials from their generating functions using 
Lagrange Inversion. This bypasses the need to use the Darboux method and 
makes it possible to obtain uniformly convergent asymptotics directly from the 
generating function. We have included this in the Appendices in Subsection l6.2l 

1.1 Some previous results 

In this section we will mention some results by other authors which we will have 
occasion to use later in this paper. One of the reasons for doing this is that 
different authors have used different labelling conventions and this will enable us 
to establish a consistent notation for use when we combine results from different 
papers with mutually incompatible conventions. We will state our results using 
the conventions in j!0| . 

Several early results in the theory of quantum walks are due to Meyer |28| . 
who considered the wavefunction as a two-component vector of amplitudes of 
the particle being at point n at time t. Let 



where the chirality of the top component is labelled RIGHT and the bottom 
LEFT. At each step the chirality of the particle evolves according to a unitary 




(1) 



4 



Hadamard transformation 



\R)^±(\R) + \L)) (2) 

\L)~±{\R)-\L)), (3) 

which is why this quantum walk is sometimes called the "Hadamard" walk. The 

particle (or "walker") then moves according to its new chirality state. Therefore, 
the particle obeys the recursion relations 

* H (n, t + 1) = + 1, t) + ^*b(" - 1, t) (4) 

*i(n, t + 1) = - i*z(n + 1, t) + i*Ji(n - 1. *)• (5) 



Meyer approached this problem from the path-integral point of view, and 
obtained expressions for the ^-functions in terms of Jacobi polynomials. The 
standard notation in pQ for Jacobi polynomials is Pn^\z) but we wish to follow 
the conventions in |29j and subsequent papers that have now become standard 
in the literature on quantum walks, and use a = n/t. We will therefore use the 
notation Jq r ' s \w). 
We find, in [S] and JU] when n^t, 

Theorem 1 (Ambainis et al. [5 , after Meyer, [28 ). 

, wo f 2 " /2_1 ^V"/Z"w9 n(0) when - t<n<0 




2 " /2 ^( t +7 l ) l /2 1) ( ) w/ien ~t<n<0 

-„/2-l}(0^+l)_ i(0) < n < t 



t — n 



(7) 



, l^flf-n,*)!- 1 (8) 
\Mn,t)\ 2 = \M2~n,t)\ 2 . (9) 

Ambainis a/, use the other sign convention, so one should interchange L 
and R (or equivalently, replace n by — n) to reflect the walk before comparing 
their results with ours. This is just a relabelling, and so their results can be 
stated as in the following theorem. We will prove the above results in the form 
below 
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Theorem 2 (Ambainis et al. |5]). When n = t (mod 2) and Jq'(w) denotes 
a Jacobi polynomial, then 



fe(M)(-l) (t -" )/2 



• (_ 1 y l +i 2 -„/2 J^-j) (0) when 0<n<t 
(_ 1 )n+i 2 n/2-i jgj^J+^ x (0) wAen - 4 < n< 

(10) 



, . , ( -2-"/ 2 - 1 X ( . 1 ' n) w „ , (0) when 0< n < t 

^t_2 j 2 «/2 ij^j^^Q) w/jen - * < n < 0. 



and 



(11) 

rP R (-n, t) = + 2, 4), (12) 

(4 - n)Vi(-», *) = (-1)"(* + n)^L(n, 4). (13) 

A few Remarks: 

1. Note that Theorem [5] differs from Theorem ^ by an external phase which 
has been dropped in previous analyses of this system; we state the sym- 
metry relations for the ^-functions rather than for their moduli-squared, 
as in 0. We will discuss this in more detail below, as some properties of 
Jacobi polynomials are required for the calculation. 

2. There is a sign error in the symmetry relations for ipR and n, 4) in 
Carteret et al. ^D] (which has been corrected in the arxiv version [lip. 
The symmetry relations will be proved in Lemmas|21 Eland equation l|77f) of 
Subsection l3.1l using some integral representations of ipii(n, t) and ipLin, t). 

3. The endpoints where n = ±t have to be handled separately, see |2Hj- For 
the starting conditions 

V>l(0,0) = 1, iM0,0) = 

the wavefunctions at the end-points (where n = ±4) are 

<MM) = (-l) t+1 2-*/ 2 t = 0,1,2,..., (14) 

Mt,t) = 0, t = 1,2,3,..., (15) 

^(-4,4)=0, 4=1,2,3,..., (16) 

Vl(-M) = (-l)'2-* /2 , 4 = 0,1,2,.... (17) 

4. The two different cases in <|11[1 and (|10l) for n > and n < can be 
combined into one case for all n satisfying — 4 < n < 4. We prove this 
later using a symmetry property of the Jacobi polynomials. Our results in 
equations (|12f> and l|13f) are refinements of the corresponding relations in 
© and |JSJ, after performing the relabelling necessary to compare results 
with different sign conventions. 
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The asymptotic behaviour for the path-integral representation has been de- 
termined in Carteret et al. ^U], starting from Theorem|3 The steepest descent 
technique was used on the standard integral representation for the Jacobi poly- 
nomial. The result was uniform exact asymptotics a in the range |a| < 1 — £, 
where s is any positive number, in terms of Airy functions. 

This technique was used earlier for Jacobi polynomials by Saff and Varga [3T] 
and by Gawronkski and Sawyer |18j ; however the connection with Airy functions 
had not been recognized as far as we know. The Airy function description 
is useful for \a\ near 1/^/2 where the asymptotic behaviour changes from an 
oscillating cosine term times i -1 / 2 (for |a| < l/-\/2) to exponentially small (for 
2- 1 /2 +e< \a\ < 1-e). 

In this paper we analyze the Hadamard walk from the Schrodinger wave- 
mechanics point of view. The earliest work on this that the authors are aware 
of is that by Nayak and Vishwanath 29 . They define 

*(0,t)=X>(n,t)e M ", (18) 

n 

(where we have used the symbol 9 for the momentum instead of the k used in 
j2Hl ) so the recursion relations above becomes 

y(9,t + l) = M e y(6,t), (19) 

where 

Me = _ e io) ■ (20) 

Thus 

(0, t) = Mfo(9, 0), *(0, 0) = (1, 0) T , (21) 

where the symbol T denotes transposition. They show that the eigenvalues 
of Me are e~ luJ ° and — e luJe where ujg is the angle in [— 7r/2,7r/2] such that 
sm(uig) = (sin#)/-\/2. They also use the other sign convention, so one should 
relabel L and R as before. Their results can be stated as in the following 
theorem. 

Theorem 3 (Nayak and Vishwanath jSH]). Let a = n/t. Then 

i _i ( 1 \n+t i ft p ie 

Mn , . l±ti^f f (l + -^L=) e-<» + »- *, (23) 
2 2.TT j_n \ v l + cos 2 9 J 

We will derive the asymptotic behaviour of the ^-functions starting from 
Theorem[3]m Section[21 We will only give the complete details for the exponen- 
tial decay range, as the calculation for the oscillatory region has already been 
done by others |29l |S] . The conventional version of the method of stationary 
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phase, as used by Nayak-Vishwanath 29 , does not work in the exponentially 
small region, as the stationary points of the phase function have left the real 
line. We will show how to extend and refine this technique so that it can be 
made to work in this situation; the modification is an application of the method 
of steepest descents. 

It is not obvious that the formulae obtained by each method are the same, 
but we will prove below in Section [3] that they are. This means that precisely 
the same asymptotic behaviour can be found using both the path-integral and 
wave-mechanics descriptions of quantum mechanics. 



2 Generalizing the method of stationary phase 

The aim of this section is to extend the method of stationary phase so that it can 
cope with stationary points that occur as complex conjugate pairs on either side 
of the real axis. We start with the integral representation of ipR in Theorem [3J 
In this representation the integration is performed along the real axis. Nayak 
and Vishwanath consider the case \a\ < 1/^/2 when there are two stationary 
points (defined below) inside the interval of integration [— 7r,7r]. When we find 
the critical points of the phase function, we obtain an equation for 9 (called k 
in at the critical points as a function of a, which is 

— ct 

COS0 Q = (24) 

v l — a 

from which the critical value of u>o, call it ujg a , can be obtained using the 
Pythagoras rule (cos 2 9 + sin 2 9 = 1) and cog = arcsin 2 ^ from |29| . How- 
ever, when |a| > l/v2 this equation no longer has any real solutions and thus 
the corresponding stationary points are no longer on the real axis, see figure |3 
Therefore the standard method of stationary phase cannot provide the exact 
asymptotics. The stationary points have "moved" off the real axis and become 
a complex conjugate pair. We shall move the contour of integration to follow 
the stationary points, whilst ensuring that the contour still goes through one of 
them. Note also that the stationary points become saddle-points on leaving the 
real line; we will return to this fact in Subsection 2.2, below. 



2.1 The function uq in the complex plane 

The key to evaluating the asymptotics for these integrals lies in the behaviour 
of the phase function ojg. We will therefore begin by describing the analytic 
properties of the function tug in the complex plane, in particular in the strip 
— 7r < $19 < 7T. We will need this information when we replace the initial 
interval of integration by a contour in the complex plane, as explained in the next 
subsection. We will also need estimates of ujg at +oo in this strip to show that 
our new contour integral converges. The singular points of tog = arcsin 
are found from the equations ^y^- — ±1. When we write 9 = u + iv, with 
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Figure 2: The behaviour of ujg as |q| > l/\/~2. As \a\ moves beyond the critical value 
of 1/ \/~2, the phase function tug becomes imaginary. The behaviour of the real part is 
shown in the solid line and that of the nascent imaginary part is shown in the dotted 
line. Multiplying this imaginary phase function by i gives a simple exponential decay. 



u G [— 7r, 7r] and v € R, we conclude from the equations sin(u + iv) — ±\/2~, 
where 

sm(u + iv) = sin u cosh v + i cos u sinh v (25) 

that u — ±^7r and coshu = y2 (or u = ±arcsinh(l)). Because of conjugation 
and symmetry there are four singular points, namely ±^7r + iarcsinh(l) and 
±i?r — iarcsinh(l). On the four half lines ±^7r + iv with \v\ > arcsinh(l) the 
function sin(u + iv) is real, and has an absolute value greater than y/2. These 
four half- lines are taken as branch cuts of the multi- valued function wq. They 
correspond with the two branch cuts of the function arcsinz = arcsin(a; + iy) in 
the z— plane from x = ±1 to x — ±oo, with y = 0. 

The strip —it < u < tt that is delineated by these branch cuts is the prin- 
cipal Riemann sheet on which tug is analytic and single- valued. We consider 
the principal branch of this function that is real on [— tt, tt] and continuously 
extended on the principal sheet. Since the function is periodic, the same holds 
for the other strips [kir, (k + 2)ir], k 6 7L. 

In FigureEHwe show the conformal mapping by u>e from the strip — 7r < u < tt. 
We show the images of a number of lines, where we concentrate on v > 0. For 
»<0a similar picture can be given. We observe the following useful facts. 

1. The image of the interval [0,7r] must go around a branch cut because 
Lug is not single valued on this interval; the image point A is given by 
A = arcsin -4= . 
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Figure 3: Conformal mapping tug = arcsin At the left is the 9— plane, 9 = 

u + iv, — n < u < 7r, with the four branch cuts. At the right is the tog — plane; only the 
image of the half-strip —n<u<iv, v > Ois shown. 



2. The points D and F are on different sides of the branch cut; the loop DEF 
around the branch cut is mapped to the vertical DEF, and the same goes 
for the three points D'E'F'. 

3. On the positive imaginary axis u = 0, v > uig has the form (cf. I|25|ll 

sinh v / sinh 2 v 



sinh v 

ujg = larcsmh ^— = I m i — 

V2 \ y/2 



+ 1 , v > 0. (26) 



4. On the half-lines u — ±7r, v > 0, ujg has the form 

sinh v / sinh 2 v 



sinhu 

ujg = — iarcsmh = —i In 



V2 



V2 



v > 0. (27) 



5. For on the vertical ED' we can write uig in the form 



I , / coshw /cosh 2 i; \ , r- 

o.v, = -7T + » In [ + W — 1 I , cosht; > V2. (28) 



6. For 8 on ED we must choose the negative square root, which gives 



I , / coshw /cosh v \ /- 
= -7T - iln I + y — 1 I , coshv > V2. (29) 



We will also need the following lemma, in order to bound ujg as the imaginary 
part of 9 tends to +00, and hence guarantee that these integrals will converge. 
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Lemma 1. If 9 = u + iv,v > 0, then, as v — ► +oo 

\0(e- Vt ), iTT < | U | < 7T. 1 ' 



e 



Proof. Since u>g = arcsin (^/?fj we have sin ujg = sin#/\/2; 



e 



.sin 



coswe — i — p-. (31) 



V2 
Hence 

/ sin 2 9 sin 6> . , 

e •= ± v 1 -— (32) 

where the ± sign in front of the first term has yet to be determined. For small 
values of 9, it is obvious that we should select the + sign, because both sides of 
equation i|32|) have to approach unity as 9 — > 0. In fact, the + sign should be 
chosen throughout the strip \u\ < hit. This is because 

_ sin^ = 3 + cos 20 _ 1 24fl 

2 4 ~ 8 [ ' 

as — ► +ioo, we conclude that in the strip |m| < ^tt 

1 -is -sinfi 1 1 n - %e 



=e — z 



e" 1 " = Q{e v ), (34) 



2^2 V2 V2 

as u — * +00. Observe that this is in agreement with the behaviour of ujg on the 
positive imaginary axis, as given in equation l|26|l . It also agrees with the confor- 
mal mapping shown in Figurc|21 where we see that the domain AEFGF' E' A' OA 
is mapped to Swg > 0. The figure also shows that the domains ABODE A and 
A' B' C D' E' A' are mapped to ^sug < 0. This corresponds to choosing the neg- 
ative values for 9 in in these domains; this gives the estimate in the second 
line of i|3U|) . This proves the lemma. □ 

Now that we have established the behaviour of the phase function u>$, we can 
proceed to choose an appropriate contour of integration. 



2.2 Saddle-point analysis 

We will now obtain an asymptotic approximation for ipR in the exponentially 
small range outside the main peaks. To find convenient locations for the contours 
of integration with respect to the stationary points, we will begin with the 
integral in ll'L'l) of Theorem [3] for iPr(2 — n,t), and use the symmetry rule for 
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ipn(n,t) (cf. Theorem|2J) to obtain the result for ipn(n,t). So, our starting point 



is 

^(n, 1) = (-1)"+Vfi(2 - n, t) = L_J / -i^-e a) t ^ 

2n J—* v 1 + cos 9 

(35) 

where a — n/t. We have dropped the factor 2 ' , because we always can 

assume that n and t have the same parity; the wavefunction is identically zero 
otherwise as the walker must always move at each time-step. 

We first locate the stationary points or saddle-points in the traditional way, 
that is, we solve the equation 

a= = ^ s9 ^ = 1 cos 9 a . (36) 
< W ! l-su?e a /2 Vl + co8^ a 

Note that in (|36|l cos^q, and a have the same sign (and a is positive). This 
gives 

cos9 a = ±-=^==. (37) 
V 1 — a 

Thus, when a < l/v2, this gives two real stationary points 

a = ± arccos f ~/= j ( 38 ) 

which are used in the stationary phase method in [2H| ■ If 1/ \/2 < a < 1 these 
points are purely imaginary, and they are given by 



ztiarccosh ( — ^=^= I . (39) 



When 1/V2 < |a| < 1 we shift the contour in the integral representation of 
the ipR given in l|35|) off the real axis to the segments shown in Figure [S] Our 
modified stationary phase method is in fact a version of the method of steepest 
descents. 

The contour of integration goes through the saddle-point on the positive 
imaginary axis, that is, through 9 a — iarccosh(a/vl — c?) and fixes the imag- 
inary part of iuig — i9a. This is equivalent to fixing the real part of we — 9a. 
We proceed as follows. Consider the integral along the contour in Figure |3J 
We can make this into a closed contour by adding in segments from — tt — > tt, 
tt — > tt + ioo and — tt + ioo — > — tt, thus obtaining an integral over 

( — 7T, 71") U (7T, 7T + ioo) U (7T + IOO, #q,) U (+#a 7 + *Oo) U ( — 7T + ZOO , —71"). (40) 



The singular points of \/l + cos 2 9 follow from solving cos 2 9 = — 1, which 
gives 9 — ±7r/2 ± iarcsinh(l) (see the open dots in Figure EJ. We avoid the 
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Figure 4: The saddle-point contour for the integral in l|M5|l for the oscillatory range 
< a < l/v2- The interval [— it, it] can be replaced by a path that runs from — it to 
— 7r + ioo; from that point through the saddle-point at the negative real axis to —ioo 
and from that point through the saddle-point at the positive real axis to +n + ioo and 
then to +7r. Note that the contributions from the vertical half lines cancel each other 
out. On the contour shown, the imaginary part of the phase function —i{iog — 9a)t 
is constant (equal to 0, in fact). The real part tends to — oo in the valleys at infinity, 
and has a maximum at the saddlepoints (black dots). 

singularities and branch cuts of the square root in the integrand and of the 
function ujq (these singularities are the same for both functions; see Figure |3J). 
The integrand is then analytic around and inside the contour in l|40|) . so the 
integral around the contour is zero. The integrals over the curves indicated 
below are therefore equal, 

(— 7r, 7r) = — (it, tt + ioo) U (tt + ioo, 6 a ) U (+0 a , — 7r + ioo) U (— 7r + ioo, — it). (41) 



The steepest-descent curves for this integral are shown in Figure [SJ Fur- 
thermore the periodicity of the integrand modulo 2ir means that the segment 
integrals from 7r to n + ioo and from — 7r + ioo to — tt cancel, so these are not 
shown in Figure El Thus, the integral from — it to it equals the integral along 
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Figure 5: The saddle-point contour for the integral in l|M50 for the exponential-decay 
range, where 1 > a > l/v2 and the countour runs from — it to — it + ioo. It continues 
from that point through the saddle-point at the positive imaginary axis to +7r + zoo and 
from that point to +7r. Again, the contributions from the vertical half-lines cancel each 
other. On the contour shown, the imaginary part of the phase function —i(u)e — 6a)t 
is constant (equal to 0, in fact). The real part tends to — oo in the valleys at infinity, 
and has a maximum at the saddle-point (black dot on the positive imaginary axis). 

the contour from — it + ioo to ir + ioo through a . From LemmaClwe conclude 
that 

as v — > +oo in the strips —tt<u<—^tt and ^ir < u < n. Hence, convergence 
at infinity on a contour as shown in Figure [S] is guaranteed. 



2.3 Evaluating the main contribution 

We evaluate e~ lUa+ aCl for the saddle-point on the positive imaginary axis, that 
is, for 

a 



9 a = zarccosh 



Now, with x = a/ \/l — a 2 , we obtain 



(43) 



iarccosh x = i \n(x + (x 2 — l) 1 ^ 2 ) = iln 



V2c 



s/T 



(44) 
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Thus 



Vl - a 2 y y ya + V2a 2 - 1 
Let us now consider uj a — arcsin ( sl ^|° ^ . We have 



cos0 a = ^==, (46) 



so therefore 



/ TT7 /l - 2a 2 /2a 2 - 1 . , 

sin d Q = y/i- a y(i - « 2 ) = y -y— = *y -y— ( 47 ) 

Thus we obtain 

^ Q = iarcsinh ( ^i-J ) ■ ( 48 ) 



Now arcsinhx = ln(x + va?+T), so 



,2 i \ 1/2 



, , v 2a 2 - 1 / 2a 2 - 1 
Ll>„ = i In — — + — ; ^r- + 1 



i In 



/ V2a 2 -1 (2a 2 - 1 + 2- 2a 2 ) 1 / 2 ' 
\yWf^ + ^/2(l-a 2 ) < 



.!„ V (49) 



V2vr 



a' 



Thus 



and hence 



exp ln V2vi-^ r V2vi-^ ■ (50) 



a/1 - a 2 \ 1 + V2a 2 - 1 



+ V2a 2 - 1 / V%VT 



a 



(51) 



Since we — 9 a a is an odd function of 9 a , we will obtain the reciprocal of this 
result for the saddle-point —i9 a . The saddle-point on the positive imaginary 
axis (when a > 1/V2) is indeed the relevant one as we will see. The exact 
saddle-point contours are shown in Figures 0] and [S] 

Now that we have established these preliminary results, we can proceed to 
prove 
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Theorem 4. Ifl/y/2 + E < \a\ < 1 - e, fften 

(a + V2a 2 - 1) t-V2 / / ^/I - ^ \" i + V2a 2 -l~ 



^2it(1- a 2 )V2a 2 - 1 



a + V2a 2 - 1 / 2-s/l - a 2 



ipL(n,t) 



(52) 

(-1)™(1 -a)*- 1 ^ // VT^ 2 \° l + V2a 2 -lV 



2^(1 -a 2 )^^ 2- ! V\ a + v2« J ^l/ 2vT 



(53) 



Proof. We prove the equation for ^r, the proof of the equation for ^ is very 
similar. Since tog — arcsin (^^fj > we have 

^ - - n ^ Sill L/ 2 = -«d - « 2 )v/2^T. (54) 
(1 + cos^ 0) ' 

Since we know that cos 2 8 a = a 2 /(l — a 2 ), we will obtain 

1 + cos 2 6» Q = 1/(1 - a 2 ), (55) 
and we already have (see equation JUJ) 



VI -a 2 



a + V2a 2 - 1 ' 
The standard formula from steepest descents tells us that 



(56) 



(_-\\n+l p -i6 a -i(uc,-8 a a)t / n 

Mn,t) - { 'I — = ^4/ 7 jg.. (57) 

2tt VI + cos 2 6» Q y t|w Q | 

The theorem then follows by using equations |5T)|l. JHSJ, G^t and □ 



Remark: The wave-mechanics calculation is conceptually much simpler than 
the path-integral analysis in Carteret et al. 10 . It is also simpler than the 
calculations of Chen and Ismail |12| . 

3 Equivalence of the two approaches 

We have now completed the calculation begun in 10 and obtained uniformly 
convergent asymptotics for the wavefunction via both methods. However, the 
functions ifir, and i/jr derived by each route did not appear to be the same. If 
they really were different, this would be very alarming, as it would imply either 
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that there is something wrong with Feynman's equivalence argument in |15| or, 
more likely, that there was something wrong with our calculation! 



The first thing to note is that the quantity raised to the power t in equa- 
tion i|52|l dominates the asymptotics of the logarithm of the functions i/>l and 
from the Schrodinger representation. Let us call it B(a); that is, 



,58) 



V2a 2 - 1 / V2Vl 



a 



These estimates agree with the asymptotics obtained using the method of 
Saff and Varga [3T| as used in [TJ]|, although this is not yet apparent. According 
to the calculation in |1U|. the corresponding quantity from the path- integral 
representation, namely B(a), should be 

(l-cO/2 



\ 1 + a I \ 1 — a 1 I 

which would seem to be a different function. 

The demonstration of the equivalence to the result obtained by Saff and 
Varga's method needs some identities, starting with 



1 + \j2a 2 - 1) = 2 (a 2 + ^2a 2 - 1) . (60) 



Combining 1 + 2a — \f2a 2 — 1 with the other quantities is rather fiddly (see 
below). To show that the two solutions l|58(l and 1)59(1 are equivalent, we will 
now employ the identity 



1 + 2a - v / 2a 2 - 1 _ 1 + y/2a 2 - 1 
T+a ~ a + V2a 2 -T 

which can easily be verified by cross-multiplication. It follows immediately from 
this identity that 



l + 2a- V2a-1\ / 1 + V2a 2 - 1 



1 + a J \a + \/2c 

and from equation l|60|l that 

I -a 2 \^ T± _( 2(1 -a 2 ) 



(62) 



a- 



2 



V2a 2 - lj V(l+V2a 2 -1) 



(63) 



V2VT^ \ 1 + V2a 2 - 1 

(64) 



1 + V2a 2 - 1 / V2VT 



a 



2 
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Thus 



1 + 2a - y/2a 2 - 1 \ / 1 - a 



\ l + a J \a 2 - V2a 2 - 1 

so that the formulas for B(a) and B(a), in lfKH| and respectively are indeed 
equivalent. 

In fact, the representations for the two functions are completely equivalent, 
but some of the steps required to prove this need some rather subtle calculations 
involving special functions, as we will now demonstrate. These technical lemmas 
will also allow us to rederive the symmetry relations in the wave-mechanics pic- 
ture that were first proved (for the path- integral picture) in [5]. For convenience 
we recall the integral representations of Theorem 

*^«»-5£( 1 + 7IT=^) (66) 

*^- t) = hL7T$m^- + " Vd6 - (67) 



where a = n/t. Note that we have omitted the factors ~ 2 since n and t 

must have the same parity because the walker must move at each time-step. 

From these integral representations obtained in the Schrodinger picture |29| . 
we prove in this section the symmetry relations for ipL and tpR and the relations 
first proved for the Jacobi polynomials. These results could previously only be 
proved in the path- integral picture [S]. 

3.1 Symmetry properties 

Lemma 2. The function tpn of 1671) satisfies the symmetry relation 

<M-n, t) = (-l) n+1 t/j R (n + 2, t). (68) 

This is the symmetry relation in Theorem |21 for ipn(n,t). 

Proof. We have, from equation (|67|l 

IT cos((n-l)fl)cosM) XT sin((n-l) fl)8i nM) 

W Vl + COS 2 6 7T J VI + cos 2 6 

(69) 

The first integral vanishes when n is even, the second one when n is odd. To 
verify this, split [0,ir] = [0, ^tt] U [^tt, it] and write on the second interval 9 = 
7r — 6'. We conclude that ipn(l — n,t) — i/jr(1 + n,t) when n is odd, and 
iPr(1 — n,t) = —i()r(1 + n, t) when n is even. This proves the lemma. □ 
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Lemma 3. The function ij)L of equation (|66[) satisfies the symmetry relation 

(t-n)M~n,t) = (-l) n (t + n)Mn,t). (70) 
This is the symmetry relation in Theorem [2 for ipi,(n,t). 
Proof. We have from (|66[l 

Mn, t) = Mn, t) + ±- f 006 g e -i(t. e+ ne) # (n) 
2n J--* V 1 + cos 2 6 



where 

(72) 

l Z" 7 ' 

cos(n0) cos(tug) d9 / sin(n^) sin(t(jfl) (73) 



27T ./_„ 



The first integral vanishes when n is odd, the second one when n is even. This 
gives the symmetry relation 

M-n,t) = {-l) n Mn,t). (74) 

Next observe that 

due _ cos 8 

d9 ~ Vl + cos 2 0' 1 ' 

and that an integration by parts in the integral in equation (|71|l gives us that 

t — n ~ 

^ L {n,t) = — —if> L (n,t). (76) 
Finally we use equations (|74|l and l|76|) to complete the proof of the lemma. □ 

Remark: The symmetry relations for ipi, and ipR also follow from the following 
property of the Jacobi polynomials: 

This formula follows from the representation of the Jacobi polynomial in terms 
of the hypergeometric function (cf. |341 p. 151, (6.35)]) in combination with a 
functional relation of this function (third line of equation (5.5) in [341 p. 110]). 
The result in (|77() combines the first case in (|10fl with the second case, and it 
also implies the symmetry rule for ipR, and similarly for (|1 1|) . 
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3.2 The ^—functions in terms of Jacobi polynomials 

We will now prove that the ip— functions with the integral representations given 
in (|66|) and i|t)7|) can be written in terms of the Jacobi polynomials as in The- 
orem |2 This step will require the use of generating functions. These provide 
a method for writing a series as the coefficients of a formal power series, in a 
dummy variable, z, for ease of manipulation. The powers of z are then the 
summation variable for the series. For a basic introduction to the theory of 
generating functions, see [SSj; for an advanced treatment, see |T5] . 

We will use these generating functions to give us an exact representation of 
the ^ functions as functions of t, thus each labelled term in the series will be the 
function for that value of t. Therefore, our approach will be based on generating 
functions that contain the ip— functions with t as the summation variable. We 
only consider sums of ^—functions with n and t having the same parity. 

3.3 Some generating functions for tjj 

For convenience, we will define the following generating functions, which can be 
obtained from the Schrodinger representation of the wavefunction. 

Theorem 5. Consider the generating functions for \z\ < 1: 

oo 

F ro (z) = ^Vfl(2m + l,2i+l> t , (78) 
*=o 

oo 

) = ^^(2m,2i)z 4 , (79) 
t=o 

oo 

H m (z) = J2$ L (2m+l,2t+l)z t , (80) 
t=o 

oo 

I m (z) = Y.^L{2m,2t)z\ (81) 

4=0 

where ipi,(n,t) is defined in (|72|l . After the summations have been performed 
these functions become, respectively, 

om— \ ~ro 

F m (z) = -== , m= 0,1,2,..., (82) 

v / TT^(i - z + VTT^) 2m 

QTti—l gin 

G m (z) = ; , m= 1,2.3..., (83) 

Gq(z) = ; (84) 
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H m {z) 



2^-5(1 + z )z 



m= 0,1,2,... (85) 



I m (z) 



v / TT^(i - z + VTT^) 2m+1 ' 

2 m - 1 (l + z)z m 



m= 1,2,3..., (86) 



Hz) 




(87) 



y/l + Z 2 Vl + Z 2 (l-Z + y/l + Z 2 ) 



These summations can be done using some fiddly, but essentially mechanical 
manipulations; we have included detailed proofs for F m {z) and a sketch of that 
for G m (z) in an appendix, in Subsection 16.11 The other generating functions 
can be obtained via similar constructions. 

For H m (z) we obtain 



It remains to compare these with some generating functions for Jacobi poly- 
nomials. For an introduction to the theory of Jacobi polynomials we refer the 
reader to |3U1 \7\ and ■ This correspondence between the two sets of generat- 
ing functions forms the central plank of the Feynman equivalence between the 
two representations of the wave-function for this system. 

3.4 Comparing the generating functions for i\) 

We now compare the generating functions l|82|) - (|87|l from the wave-mechanics 
representation, with the generating function of the Jacobi polynomials (cf. |7| 
p. 298]) from the path-integral representation to complete our proof of the 
equivalence of the two approaches: 




(88) 



(89) 



(90) 





(92) 



where R = \/l — 2xz + z 2 , which for x = becomes 

°° , . cyr+s 




(93) 
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By applying the Cauchy integral formula we find that 

j(^) (0) = ^_i g g (94) 

where the integral is taken over a circle with radius less than unity. 
Lemma 4. For n = we have 

W0,0) = 0, (95) 
^(0,2t) = ±j£f(0) - ^(-ir'jfi^O), t = 2,4,6, ... . (96) 
For < 7i < t, n and t having the same parity, we have 

i, R {n,t) = 2-*(_l)*i a (_l)»+ij2£- 1) (0). (97) 

2 

This gives the first case of equation (|10[) from the wave-mechanics representa- 
tion, which was originally obtained via the path-integral method. 

Proof. From ffltft and (JH2J we find, as in IpEI) . 

ib R (2m + l,2t + 1) = <b —== ; — r . (98) 

We now compare (1941 with l|98[) and take u = 2m, u = 0, and s — t ~ m. This 
gives 

i> R (2m + 1, 2i + 1) = 2- m -* J t ( !™' 0) (0), < m < t. (99) 
Using the symmetry rule for the Jacobi polynomials 

j<r\-x) = (-i)«j<rHx), (ioo) 

which follows from 1)92(1 by putting x — > — x, z — > — z, we find 

V'fl(2m+l,2y;+l) = 2- m -^(-l) t - m J t ( ^ l m) (0), < m < t. (101) 

For the even case, we obtain from (|83[) . and 11 QUI) 

V^(2m, 2t) = -2- m J t ( !™- ll0) (0) - -2- m (-l)*- m J^ m_1) (0), < m < f, 

(102) 

and from J7SJ and we obtain lj9"5)l . Combining (|101fl and l|102|l in one 

formula gives l|97|) . This proves the lemma. □ 

Lemma 5. For n = we have 

1^(0,0) = 1, (103) 
tf L (0, i) = J t % 0) (0) - \j^\{<S), t = 2,4, 6, . . . . (104) 

For < n < t, n and t having the same parity, we have 

Mn,t) = 2-«/ 2 -i ( _ 1) ( t -)/ 2 -i J a-) )/2 _ i(0 ). (105) 
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This gives the first case of from the path-integral representation, now ob- 
tained via wave-mechanics. 



Proof. We obtain from (JSOJ) and JSSJl 

M2m + l,2t+l) = 2" m -i [ (0) + (0)] (106) 

£ + TO + 1 ' ™ W ' V ' 

where we have used the relation for the Jacobi polynomials (cf. 1;, p. 782, 
(22.7.19)]) 

(u + v + 2k)4 u ' v ~ 1) (x) = {u + v + k)4 u ' v) (x) + (u + k)j£*\x). (108) 
We use ff7)l. ifTBl and (JTHOl, and obtain 

^ L (2TO + l,2t+l) = 2- m -i(-l) t - m - 1 J t ( ^™l 1) (0). (109) 
For the even case we obtain from (|81|) and H86|) 

M2m, 2t) = 2— 1 [j t ( !™' 0) (0) + 4^(0)1 (110) 



= 2-™- 1 (-l)*-™- 1 ^-J t ( ^™ ) 1 (0), (111) 
where we used (|108f) . (I77jl and j 100(1 . By using JZHJ we obtain 

^(2to, 2t) = 2 -m-i(_ 1)t -m-i jJW (0)) m=l,2,3 (112) 

From (|HU and JgZJ we obtain (|103|) . Combining (|109|l and (|112|l into a single 
formula gives equation (|105fl . This proves the lemma. □ 

As we have now established that both sets of generating functions match, this 
completes the proof of the equivalence of the results obtained via the path- 
integral and wave- mechanics representations. It should be noted that we have 
proved the two representations are exactly equivalent for all time, as opposed 
to being only asymptotically equivalent in the long-time limit. It is one of the 
curious features of generating function methods that they can be used to prove 
the existence of a one-to-one correspondence between the two sets (i.e., our 
V'-functions) counted by the two series, without actually finding the explicit 
bijection. 

3.5 Summary of the equivalence results 

While the coined quantum walk can be thought of as a quantum analogue of the 
discrete-time classical random walk it should be noted that the quantum 

model inherits its discrete time parameter directly from the classical model; the 
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discreteness was not introduced by hand as part of the quantization procedure. 
Also, we have not defined a Hamiltonian for this system at all, so the problem 
of ambiguities in the time derivatives of the action does not arise. We have now 
shown the full Feynman equivalence for this system, though some results seem 
easier to derive in one approach than in the other. 

• We have obtained the symmetry rules directly from the integral represen- 
tations for the ^-functions. It is not necessary to represent the ^-functions 
as Jacobi poynomials and then use the symmetry properties of Jacobi 
polynomials (as was done in to obtain this result. 

• The relations between the ^-functions and the Jacobi polynomials have 
been obtained directly from the integral representations, though we needed 
to develop some technical tools for this. These consisted of a few extra 
properties of the Jacobi polynomials, and the generating functions con- 
taining the ^-functions. The proofs using these generating functions are 
conceptually straightforward, although the details are quite technical. The 
Feynman path-integral approach of [Ej (which is a finite sum here) would 
seem to be simpler if these relations are all one wants. 

• We were able to establish some new expressions for the values of the two 
components of the ^-function at n = as a function of time, in equations 
(|95l) and l|103|l . which were not known before. 

4 Physical interpretation of these results 

Now that we have established that the rather counter-intuitive results obtained 
from the wave-mechanics picture really are equivalent to the Airy functions ob- 
tained from the path-integral approach, we are left with the little mystery of 
their physical interpretation. In this region of exponential decay these waves 
have complex wavenumbers. This phenomenon is called evanescence. And 
herein lies the mystery; the conventional wisdom is that evanescent waves are 
only ever seen in the presence of absorbing media, such as light waves being 
absorbed into a conducting surface, but there is no such surface here and the 
evolution is unitary, by the initial assumptions that went into constructing the 
model. In fact, the phenomenon of evanescence is rather more widespread; it 
occurs in a great many systems if you know where to look. In a pioneering 
paper in the early 1990s, Michael Berry showed that evanescent behaviour is 
much more common than had been previously thought, after being inspired by 
some work by Aharonov, Anandan, Popescu and Vaidman in [3]. Berry gave a 
detailed discussion of how this phenomenon occurs in optics, at the edges of the 
almost ubiquitous "Gaussian" beams in [5]. 

The wavefunction for this system tends to an Airy function in the asymptotic 
limit, as was proved analytically in 10 . We evaluated the integrals in the 
path-integral picture using the method of steepest descents, which in this case 
featured a pair of coalescing saddle-points |101 I27| . Since then, various authors 
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have discussed the connection between the discrete walk on the infinite line 
and interference phenomena in the quantum optics of dispersive media. This 
connection was first described by Knight, Roldan and Sipe in a series of papers 
|25l 1241 l2l)| and further clarified by Kendon and Sanders in [23| ■ 

As we enter the exponential decay region, the two original stationary points 
of the phase function merge and then two new saddle-points are born, which 
move off the real axis as a complex conjugate pair. The behaviour of the mo- 
mentum closely follows that of uig, which was plotted in Figure |2] In the wave- 
mechanics picture, we found that the momentum becomes purely imaginary in 
the exponential decay region; indeed, the techniques we used to evaluate the 
integral relied on this fact. So, the behaviour of the walk in the exponential 
decay region is a pure exponential decay; there is no oscillatory behaviour. 

Within the interpretation begun by Knight et al, it was first suggested to 
us by Achim Kempf 221 that the specific evanescence phenomena that we have 
discussed in this paper are analogues of what are called the Sommerfeld and 
Brillouin precursors. Specifically, the exponential decay region can be identified 
with the Sommerfeld precursor (see for example, |17|1 and the distinctive 
peaks in the probability distribution would be an example of the Brillouin 
Precursor (see for example ^H])- Our results here and the previous results in 
[TU| provide the first analytic evidence for this identification. 

5 Conclusions 

In this paper we have completed the analysis begun in JU|, thus meeting the 
challenge made in |5| to prove all their theorems about the unrestricted quantum 
walk on the line in both the path-integral and wave-mechanics representations. 
We have also proved some additional identities that we believe to be novel. In 
the course of doing this, we have had to generalise the method of stationary 
phase in a way that may have applications beyond this problem. We have also 
proved the exact Feynman equivalence between the two representations directly, 
by reducing the problem to purely combinatorical constructions which may also 
be of wider interest. 

Lastly, we have supplied a physical interpretation for our results, in terms 
of certain evanescent phenomena from the quantum optics of dispersive media. 
This interpretation is a somewhat counter-intuitive one, as it would seem to 
require an effective dissipation that acts on the walker in a way that is analogous 
to the effect of a dielectric medium on light, despite the fact that the evolution 
of the system is unitary by assumption. 
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6 Appendices 

6.1 Construction of the generating functions 

Here we give the details of the construction for the generating functions from 
Subsection l3.3l 
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6.1.1 Proof of the construction for F m (z) 

Proof. We give a detailed proof for F m (z). We substitute equation (|rJ7|l into 
equation (|78|l and obtain 



7T -i(2m9+u ( ) -i 

, ;— - dff 

VI + cos 2 9 1 - ze~ 2lu e 

l nn e - t (2m8+ U9 ) 1 _ 2 C0S 2u)g ~ IZ sin 2u>g 



2tt J_ n VI + cos 2 6> 1 - 2z cos 2w + z 2 
1 {cos(2?ti6' + uig) — i sm(2m8 + ujg)}{l — z cos 2uig — iz sin2cjg} 



d0 

in 0/,in\ 

d9. 

27T Vl + cos 2 (1 - 2z cos 2u>g + z 2 ) 

(113) 

Since tog is an odd function, the imaginary parts will vanish. The interval [— w, tt] 
can be reduced to [0, w] and we obtain 

. . 1 cos(2m0 + ujg)(l — z cos 2uig) — z sin2wg sin(2m0 + uig) 

r m U = — / ; = do 

1 Vl + cos 2 0(l-2zcos2w e + z 2 ) 



7T 



(114) 



COs(2m9 + Uig) 

do 



tt Jo Vl + cos 2 (1 - 2z cos 2ug + z 2 ) 

z cos(2to0 + tug) cos(2uig) + sin(2m0 + uig) sin 2wg 

/ = ctw. (115) 

ttJo V 1 + cos 2 (0) (1 - 2zcos2cj 6 + z 2 ) 

The numerator of the second integral can be written as cos(2m0 — u)g), and it 
follows that 



F m ( Z ) = - r 

IT JO 



cos(2m0 + ujg) 

do 



Vl + cos 2 (1 - 2z cos 2uig + z 2 ) 

_z_ r co S (2m9-u e ) 

Wo Vl + cos 2 9 (1 - 2z cos 2w + z 2 ) 

We can now use simple trigonometric identities for the cosines in the numerators, 
to obtain 

If* cos 2m9 cos uig — sin 2m9 sin uig 

F m (z) — — I d9 

^ Jo Vl + cos 2 6 (1 - 2z cos 2ujg + z 2 ) 

z f* cos 2m9 cos tog + sin 2m9 sin ajg 

/ = a0. (117) 

Wo Vl + cos 2 9 (1 - 2z cos 2w 9 + z 2 ) 

Now we use sin uig — (sin 0) / V2, and observe that the terms with the sine 
functions do not contribute to the integrals; this follows easily by performing 

the transformation 9 = 9' + Using also coswg = + cos2 we obtain 

, s l-z r cos(2m0) ,„ , s 

F m (z) = ■= I — -d9. 118 

V ' nV2 Jo l-2zcos2w e + z 2 1 ; 
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For the final step we use formula (3.615) (1) of Gradshteyn and Ryzhik |20|. 
that is, 



2 in 



^ cos (2m6) , n (-l) m 7r (1 - VI - a 2 ) , , 

i f-dfl = - — , a 2 < 1, m = 0,1,2,... 

1- a 2 sin 2 6- 2VT^ a 2m 'II' 

(H9) 

We observe that cos2wg = cos 2 0, and take a 2 = —2z/(l — z) 2 . This gives the 
expression in (|52l . as advertised. 

6.1.2 Outline of the proof for G m (z) 

The proof of equation (|83|l for G m (z) uses essentially the same manipulations: 



e i(l-2m)S 2 

27r J_ 7r Vl + cos 2 6» l-ze 

1 cos(2m — 1)0(1 — zcos2oj£/) — zsin(2m — 1)6* sin2w, 



d6>. 

27r J_ 7r VI + cos 2 6» (1 - 2zcos2cj e + z 2 ) 



(121) 



The second integral can be broken into two terms as follows: 



^ / x 1 r cos(2m-l)0 ,„ 

G m (z) = - / '- d9 122 

Wo Vl + cos 2 6l(l-2zcos2w ( , + z 2 ) 

z /"^ cos(2m — 1)6* cos 2w6i + sin(2m — 1)6* sin 2oje 

/ — da. (123) 

Wo \/l + cos 2 6> (1 - 2z cos 2w e + z 2 ) 

The first integral vanishes (substitute 9 = 9' + hn). The same holds for the 
contributions from the cosine terms in the second integral. This gives 



z 



r sin(2m-l)0sin0 

G m z = / — ^ ' dd. 124 

7T 7n 1 — 2z COS 2u)g + z z 



Using the fact that 



cos2m6> = cos(2m - 1)6* cos 9 + sin(2m - 1)9 sin 9 (125) 

cos(2m - 2)9 = cos(2m - 1)6* cos 9 - sin(2m - 1)9 sin 9 (126) 

we can then write 

„ z r cos 2m0-cos(2m- 2)6> M . 
27r ,y 1 — 2zcos2we + ^ 

and we see (cf. Q118[l l that we can express G m (z) in terms of two F m (z) func- 
tions. This gives equation Q83[l. □ 
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6.2 Lagrange inversion asymptotics 



Suppose we have an unknown function w that we assume can be written as an 
(as yet unknown) power series. All we know about this power series is that it 
can be written as a recursion relation 

to = z(p(w), ip(0) f 0, (128) 

where tp(w) is some generating function that is defined as an implicit function 
of z. We would like to find w as an explicit function of z, so we can express 
some other function f(w) as an explicit power series in z. Lagrange Inversion 
enables us to do this, and tells us that we can write 

/H = 2-[A"- 1 ]/(A)^(A), (129) 

n>l 

where A is a dummy variable, ' denotes differentiation with respect to A and the 
square brackets [x n ] is the "Goulden- Jackson" notation ^5] for the coefficient 
of the term in x n . 



Here are two very simple examples to illustrate the use of Lagrange Inversion. 
Suppose we were given w defined only to be a solution of the equation 

to = ze w , (130) 

and we know nothing else about it. Suppose we just want to obtain a power 
series for /(to) = to, where <p(w) = e w . Then the formula above becomes 

. ill r, n ~ 1 

»=E^ = EV- (131) 

7»>1 n>l 

It follows from Stirling's formula for the factorial that this series converges for 
\z\ < 1/e. A slightly less simple example occurs if we are interested in the 
function g(w) — w 2 . Then equation 11311) would become 

j-n fn ^n— 3 

w 2 = J2 — [A n ~ 1 ]2Ae" A = 2J2 -[A"- 2 ]e" A =2^t" (132) 

n>l n>l n>2 v ' 

6.2.1 Lagrange inversion for Jacobi polynomials 

The formula of Lagrange most useful for Jacobi polynomials is from 



' /,A1 ( 133 ) 



1 — zip'(X) n! dx 

r v ; n=0 



where <p was defined in equation p28|l and x is another dummy variable. 

The asymptotics of the Jacobi polynomials have been discussed previously 
by Chen and Ismail and Ambainis et al. used those results to derive some 
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results on the asymptotics of the ^-functions. It should be noted that the 
Chen-Ismail results used the method of Darboux and are not uniform over the 
full range of a. For the rest of this section we will briefly discuss Lagrange 
inversion and show how to derive the integral representations that Carteret et 
al. used for the ^-function to derive uniformly convergent asymptotics. Chen 
and Ismail's work on the Jacobi polynomials uses the generating function 
for jjT+ a J'/5+^)(0) of Srivastava and Singhal &\ which in our notation becomes 

oo 

^jj7+<^+W (0 y = (i +u )^+ 1 (l+v)' 3 + 1 [l-au-6v-(l+a+6)uv]- 1 , (134) 

3=0 

where u = — v and u is implicitly defined as a function of z by 

-u= ^(l-u) 1+a (l + u) 1+b , (135) 

In order to obtain the asymptotics, we will need to interpret equation (|134|l in 
terms of equation l|133fl . again following Srivastava and Singhal |32) . In the case 
of i/jl, we should let (see [TU] ) 

a -7-0, P=T^> b=- . (136) 

1 — a 1 — a 

Then we should define 

^(AH^rr^l + A) 2 ^ 1 -") (137) 

and / should be defined by equation (|133|l . which we must set equal to 

/(A) = (l + A) (1+a)/(1 - a) , (138) 

following the method originally developed in [22] ■ We omit the details, since we 
only wish to use the result. Then 

4°^ / l 1 (0) = Jr°- m ^ +! ^ ) (0), (139) 

and so 

m= (l-a)t/2-l, (140) 

as required. 

We would like to obtain an integral representation of the coefficients in equa- 
tion (1 1 34(1 . If /(A) and (p(X) are analytic then this can be done using the Cauchy 
integral formula, thus 

[A"]/(A)0"(A) = -L / /(A)0"(A)A-"- 1 dX (141) 
J c 
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where C is a sufficiently small contour around the origin. Equation (|141|l is an 
example of a Rodrigues formula for a set of orthogonal polynomials. Another 
example of these appears in the method for generating orthogonal polynomials 
using Gram-Schmidt orthogonalization [jjj. So we obtain the integral represen- 
tation for the Jacobi polynomials 

jCo.Wtt-oO+flfo) = _L J (1 + A)^"^ 1 -) + A f"/(i- a )j" dA; 

(142) 

as used in [311 IT%| . This is the contour integral that Saff and Varga estimated 
using steepest descents. This is discussed in complete detail in ^H] so there is 
no need to say more here about this particular example. This example shows 
however that expressing a coefficient obtained using Lagrange inversion as a 
contour integral in this way and using steepest descents may lead to uniform 
asymptotic expansions over a wide domain. See the book [Jj by Andrews et al. 
for examples of Lagrange inversion arising in special functions. 
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